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1. Quantum Mechanics Isn't Intuitive ! 

It is apparent to anyone who thinks about it that, to a large degree, the basic concepts of 
Newtonian physics are quite intuitive, but quantum mechanics is not. My purpose in this talk is to 
introduce you to a new, much more intuitive way to understand how quantum mechanics works. 

I begin with an incredibly easy way to derive the time evolution of a Gaussian wave-packet for 
the case free and harmonic motion without any need to know the eigenstates of the Hamiltonian. 
This discussion is completely analytic and I will later use it to relate the solution for the behavior 
of the Gaussian packet to the Feynman path-integral and stationary phase approximation. It will 
be clear that using the information about the evolution of the Gaussian in this way goes far beyond 
what the stationary phase approximation tells us. 

Next, I introduce the concept of the bucket brigade approach to dealing with problems that 
cannot be handled totally analytically. This approach combines the intuition obtained in the initial 
discussion, as well as the intuition obtained from the path-integral, with simple numerical tools. 
My goal is to show that, for any specific process, there is a simple Hilbert space interpretation of 
the stationary phase approximation. I will then argue that, from the point of view of numerical ap- 
proximations, the trajectory obtained from my generalization of the stationary phase approximation 
specifies that subspace of the full Hilbert space that is needed to compute the time evolution of the 
particular state under the full Hamiltonian. The prescription I will give is totally non-perturbative 
and we will see, by the grace of Maple animations computed for the case of the anharmonic os- 
cillator Hamiltonian, that this approach allows surprisingly accurate computations to be performed 
with very little work. To view these animations go to http : \\slac.stanford.edu\ ~ niv. I think of 
this approach to the path-integral as defining what I call a guided numerical approximation scheme. 

After the discussion of the anharmonic oscillator I will turn to tunneling problems and show 
that the instanton can also be though of in the same way. I will do this for the classic problem 
of a double well potential in the extreme limit when the splitting between the two lowest levels is 
extremely small and the tunneling rate from one well to another is also very small. 

2. Gaussian Wavefunctions and the Path Integral 

Discussions of one-dimensional quantum mechanics usually begin by considering the space 
of square integrable functions \ir{x) on the interval —0° < x < 0°. This space of functions is acted 
upon two operators, p and x, where the action of these operators on a state ^f{x) is defined to be 

X = (2.1) 

pi//(x) = -^(x); (2.2) 
I dx 

i.e., X is simply multiplication by the variable x and p is differentiation with respect to the variable 
X. Given these definitions it is simple to show that x and p satisfy the commutation relation 

[x,p]=j. (2.3) 

In what follows we will devote a great deal of attention to the properties of Gaussian wave- 
packets. For our purpose a Gaussian packet of width I /y^is defined as the solution to the equation 

(jp-Fyx)vA(x) =0. (2.4) 
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To show this is the same as the usual condition is simple, but worth doing expUcitly since we will 
use it over and over again. The steps are 

{^{x) + y^if{x))=0 (2.5) 
i/^W =Ce"2^ (2.6) 
where constant C is determined by the normalization condition 

/oo 
dxY{x)w{x) = 1, (2.7) 
-oo 

that means that C = (|) 

From this point on we will write a Gaussian satisfying this equation by the symbol |0^). 

3. Shifted Gaussians 

In what follows it will be important for us to consider Gaussian wave-packets centered about 
points X 7^ 0; i.e., Gaussian packets shifted away from the origin. As is customary, these packets 
will be constructed by applying the operator 

[/(f) = ^-p-- = £iz!)!p«^~« (3.1) 

(=0 

to the state |0y) to obtain the state \xy); where x is an arbitrary number and p is the momentum 
operator. 

To see that this operation does what we want observe that 

e-m^^v^ = x-/f[p,x] + ^x2[p,[p,x]] + ... (3.2) 
= x — x (3.3) 

where all but the first two terms of the expansion vanish since the commutator of p with x is 
proportional to the unit operator. Thus, 

e-P*(/p + yx)|0^) = 0, (3.4) 
g-'P^ (/p + 7x) e'P^~ {e-'^' \0y)) = 0, (3.5) 
g-'P^" (/p + yx) e'^-' \xy) = 0. (3.6) 

Now observe that p commutes with U (x) so that Eq.^ becomes 

{ip + Y{x-x))\xy) =0. (3.7) 
Using the previous argument, thinking of \xy) as a function of x, we have 



e-i^'-'> . (3.8) 

71/ 

This, of course, is what we wished to show. 

I have only included this elementary discussion to show how one can manipulate the equation 
that defines a Gaussian packet in order to obtain a useful result. In the next section I will use the 
same sort of argument to derive the evolution in time of an arbitrary Gaussian packet. 
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4. Evolution of a Gaussian With Free Hamiltonian 

The operator form of the time-dependent Schrodinger equation, given a Hamiltonian H, says 
that the state \ '^f{t)) satisfies the equation 

i.\^f{t)) = -iuw{t)). (4.1) 

As is well known, the solution to this equation is 

IVA(O) I V^(0). (4.2) 

To derive the time evolution of a Gaussian packet we multiply Eq. 2.4 by e^"^, where H is defined 



to be 

H = ^, (4.3) 
2m 

for some mass m; i.e., we consider 

e-"^^ (ip + 7x) |0^) = (4.4) 

2 2 

e^"^' (ip + 7x) e"^" \Oy{t)) = (4.5) 

(/p(0 + 7x(0) |0r(0) = 0; (4.6) 

where the time dependent operators x(f) and p(f) are 



x(t) = e^'' ^- xe''^' =x-t-, (4.7) 

m 



p(f) = e'" ^- pe"^- =p. (4.8) 
These results follow from the definition of the exponential and the commutation relations. Substi- 



tuting this into Eq. 4.6 we obtain 



(/p + 7(x--p)|0,(0) =0 (4.9) 
m 

1 + ^^ /p + 7x^ |0^(0) = (4.10) 
(/p + 7(0x) |0^(0) = 0. (4.11) 

This, as we have already shown, means 

|0^(^))=C(0^?"^rW^^ (4.12) 

where 



7(0 = v-^- (4-13) 



The fact that 

C(0=(i , (4.14) 



Try + 

m 
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follows directly from the equation 

|(C(0.-^^«-^') = -'-C(0^^-^^^''^), (4.15) 

which is simply a differential equation for ln(C(f)). I have bothered to include the entire derivation 
of the time-dependent wave-function to show how powerful manipulating the defining equation for 
the Gaussian packet can be. Because the time dependent packet drops off as e^^-^^'^^^ where y{t) 
is complex, we will refer to this as a generalized Gaussian packet. An animation showing how this 
wave-function looks as it evolves in time is shown at the beginning of the html file that I referred 
to previously this paper. It is followed by animations showing how different coherent states evolve 
using the free Hamiltonian. 

5. The Harmonic Oscillator 

Now that we have the exact solution of the time evolution of an arbitrary Gaussian packet under 
the free Hamiltonian, let us consider the next simplest case, the harmonic oscillator Hamiltonian; 
i.e. 

// = i^ + -— x2. (5.1) 

2m 2 

Once again if we start by multiplying the defining equation of a Gaussian packet by the exponential 
of H we find 

e-'-'ff (/p + yx) |0t,) = 0, (5.2) 

or 

(/p(0 + 7x(0)|0r(0)=0, (5.3) 
where the time-dependent operators p(f) and x{t) are defined to be 

V{t)=e-''^ Ve''^ and x{t) = e-''^ x e''^ . (5.4) 

It follows immediately from this equation and the commutation relations of p and x that 

dx \ , , dp 7 . . 

— = -p? and -^ = -(o^x(t). (5.5) 
dt m dt 

Since this is simply a first order differential equation with the boundary conditions p(? = 0) = p 

and x{t = 0) = X, it has the unique solution 

x(?) = cos(g)?) x + — !— sin(«f)p (5.6) 
p(f) = cos(ft)?) p-mft)sin(ft)f)x (5.7) 



Substituting these into Eq. 5.3 we obtain, as before, 




/ / cos(a)0 + ^ sin(a)? 

h'P + 7 , , . , : ) X ) \Qy{t)) = (5.8) 

Thus we see that the evolution of a Gaussian packet in a harmonic oscillator potential is once again 
a generalized Gaussian with a y{t) that is periodic in time. Note, if and only if 7= mco, the packet 
doesn't change in time and it is simply multiplied by a phase g^"®/^; j e. when y=mco the packet 
is an eigenstate of the harmonic oscillator Hamiltonian. 
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6. Coherent States 

Previously we discussed Gaussians shifted to a mean position x. Now we will generalize the 
shifted state to one that has a non- vanishing expectation value for both x and p; i.e., consider the 
time evolution of the state 

^-/px+,px (6.1) 

This is called a coherent state. 

Clearly, using the previous arguments, multiplying this state to the left by e^"^ we obtain 

e-iv{t)x+ipAt) \0^{t)). (6.2) 
Collecting terms this can be rewritten as 

g-'PW(0+'Pda..(Ox \Q^(t)). (6.3) 

where |Oy(f)) was calculated in the preceding section and fciass(0 and yJciassCO are the solutions 
to the classical equations of motion for a particle moving in a harmonic potential that initially is 
located at the position x with momentum p. It is not an accident that this is the same trajectory 
one would obtain by doing the stationary phase approximation to the path integral. To make the 
analogy between the formula for the propagation of an arbitrary shifted Gaussian and the stationary 
phase approximation more striking let us rewrite Eq. |63| as 

g!Pclass(0 -fetes (f) g-!P-Vctes(0 + '/Ws(f)(x--Vctes(0) \Oy[t)) . (6.4) 

This form of the generalized Gaussian packet shows that the packet center moves along the clas- 
sical trajectory for a particle starting with the given initial mean position and mean momentum. 
Furthermore, it shows that, if we write the position dependent phase factor that gives the shifted 
packet the correct mean classical momentum, so that it has the value one at the packet center, then 
the entire packet is multiplied by a time dependent phase factor that is the exponential of the clas- 
sical action. These results are also what is seen in derivations of the path-integral using coherent 
states with a single fixed value for 7. What is not captured in the coherent state derivation of the 
path-integral is the fact that 7 changes in time and, in fact, becomes complex. This is why that 
approach is less powerful than what I will do next. 

7. The Bucket-Brigade Approach to the Path-Integral 

We now understand how a generahzed Gaussian packet propagates with the free or harmonic 
oscillator Hamiltonian. Now, let us spend a few moments connecting this knowledge to the usual 
derivation of the path integral using Gaussian coherent states. Most derivations begin by rewriting 
the time evolution operator as a product of the evolution operator for many small time steps and 
inserting a complete set of states between each term in the product;i.e., 

(7finalk-''^|7mit) = (7final|...|7;+l)(7;+lk-''"/"l7))(7^^ (V-D 

In Feynman's derivation of the path integral this complete set of states aie 5-functions of x or p, 
as appropriate. Later derivations used coherent states, since the shifted coherent states form an 
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over-complete basis in terms of which one can construct a resolution of the identity operator. In 
either case, after deriving this identity one customarily makes the stationary phase approximation, 
which in effect selects a single intermediate state at each step. In all of these cases, however, even 
for the exactly solvable cases of the free particle, or the harmonic oscillator, the states selected by 
this step are not good approximations to the true evolution. We have already seen that the correct 
evolution of a free particle or a particle in a harmonic oscillator potential is a generalized shifted 
Gaussian with a complex Y{t). Obviously, if we insert these states as intermediate states then the 
stationary phase approximation would produce the exact answer This observation almost brings us 
to the formulation of the bucket-brigade approach to dealing with the Schrodinger equation. The 
missing step is the observation that as the number of steps in the decomposition of the evolution 
operator increases the number of states selected by the stationary phase approximation increases 
too. However, since these states are not orthogonal to one another, the number of significantly 
linearly independent states doesn't grow in the same way. To be more precise, I define the notion 
of significantly linearly independent states as follows: let the integers M and A^, with M <N define 
two decompositions of the time interval in the decomposition. Let |]^) and \\ffj) be the two sets 
of states defined by the corresponding stationary phase condition, i.e. let them be the generalized 
Gaussian packets obtained by exactly propagating the initial state a time t/Movt/N. The larger set 
of states will not be significantly linearly independent of the smaller set if all of the larger states can 
be represented to some pre-defined accuracy as a linear combination of the smaller set of states. 

Since dividing the time interval over which the evolution is occurring into ever smaller slices 
does not lead to increasing numbers of significantly linearly independent states, it follows that one 
can describe the continuous time evolution of the initial state to arbittary accuracy be restricting 
attention to a finite dimensional sub-space of the full Hilbert space. In this section I will show that 
this is the case for free and harmonic evolution. To be precise, I will show that in order to compute 
the states £"''^1%) for all values of t between some ^initial and ?fi„ai to high accuracy it suffices to, 
given some discrete set of states \^n)^ compute the truncated operators 



and then exponentiate the finite matrix Hnm after transforming it to the orthonormal basis defined 
by the |v';,)'s. In the accompanying .html I show how well this works for various initial states 
evolving either under the free Hamiltonian or the harmonic oscillator Hamiltonian. 

8. A Non-Trivial Example: The Anharmonic Oscillator 

I have argued that the bucket-brigade idea says that, with no significant loss of accuracy, we 
can restrict attention to a relatively small subspace of Hilbert space to compute the continuous time 
evolution of a given packet, I will now show that the same is true for Hamiltonians for which the 
time evolution cannot be exactly computed. To show how this works let us begin by considering 
the case of the anharmonic oscillator; i.e., the system defined by the Hamiltonian 



and A^„m = (rn|N|ym), 



(7.2) 




(8.1) 



We will start, as before, with a Gaussian packet defined by the equation 

(jp + yx) |0^)=0. 



(8.2) 
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Then, it follows that the time-evolved packet satisfies the equation 

e-''H(/p + 7x)|T^,)=0, (8.3) 

or 

(/p(0 + 7x(0)l7(0)=0, (8.4) 

where x(f) and p(f) are defined as in Eq. |4.8| . It not possible, however, to compute x(f) or p(f) 
exactly. We can however non-perturbatively approximate the evolution of a state 

|..,-,A-,7(0>=^""''P+''"''|0r(0> (8-5) 
by evolving it with the effective quadratic Hamiltonian 

Heff(^,-,A-,7(0) = ^P;^^'^' +A(/?g(y(O)|V(x + x,0|/?g(y(O)) (8.6) 
2m 

+ X{Re{y{t))\^y{iL + Xi)\Re{y{t))){x-Xi) (8.7) 

+ X{Re{y{t))\-^\{x + Xi)\Re{y{t))){x-Xif. (8.8) 

Note that in order to guarantee that the Hamiltonian is hermitian the expectation values are com- 
puted for Gaussian packets where y{t) is replaced by a Gaussian with the same where y{t) is 
replaced by the real part of y{t). Thus, with this in mind we define the iterative procedure where 
we begin with a generahzed Gaussian 

I \lfn) = Cn e'P" ^'-'"^ e--2 ^" (8.9) 

and evolve it with a quadratic Hamiltonian with the generic form 

Heff= hVn-Fn{x-Xn) + —r—. (8.10) 

2m 2 
Applying the formulas we already derived for the evolution of a generalized Gaussian in a harmonic 
potential we see that we get a new generalized Gaussian of the form 

\Y„+i) = C„+ie'^^^'^"+' {x-{x„+x{St)„+i]) ^-^^ (x-{x„+x(St)„+i))- (8.11) 

where the quantities appearing in this equation are given by the formulas 



1 



and 



C0&{(0n5t) + sin{(Ont) 



x(dt)n+i = ^^(\-cos( (On 5t)) + -^sm((0ndt) (8.13) 

f 

p(5t)n+i = cos(g)„ 5t) p„ + — sin(G)„ 5t) (8.14) 

(On 

cos(a)„ 50 + sin(a)„ 5t) 
cos(w„ 5t) + j^sm{(On5t) 
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Note, that at the n step the Hamiltonian parameters, for the case of the anharmonic oscillator are 
defined by the equations 



One repeats this process again and again, recomputing the effective Hamiltonian at each stage, 
and obtains a basis for the truncated Hilbert space. To compute the continuous time evolution 
of the initial state we exponentiate the truncated Hamiltonian, obtained by computing the matrix 
elements of the exact Hamiltonian between these Gaussian packets. Since these are Gaussians, 
this is easy to do. The comparison of this approximate computation and the exact result obtained 
from a brute force numerical approximation is shown in the previously referred to html file. That 
file contains many animations showing how a free Gaussian packet evolves in time under various 
circumstances, as well as how a Gaussian in a harmonic potential evolves in time. It also contains 
animations that compare the evolution using discretized, bucket-brigade states to exact solutions, 
both for the exactly solvable systems, as well as for the anharmonic oscillator. You will see that, in 
all cases, the agreement between the real and imaginary part of the wave-fucntion for approximate 
and exact calculation is quite remarkable. 

9. Tlinneling and Instantons 

The final issue I want to touch upon is tunneling, which is important both for problems related 
to tunneling between different minima of a potential, and to general problems of scattering from a 
non-square barrier. 

Consider the Hamiltonian for a particle in a double- well potential 



If we now attempt to find stationary Gaussians there will be two solutions, one in each well. As 
before the classical moments, p of each solution must be set to zero. Next, as is shown in the .html 
file, the condition that the force, the expectation value of the derivative of the potential, vanish 
means that in the left-hand well the Gaussian is shifted sUghtly to the right of the minimum and 
in the right-hand well the Gaussian is shifted slightly to the left. The parameter 7 for each of the 
Gaussian's doesn't evolve in time determines 7 in terms of the appropriate expectation value of the 
second derivative of the potential. If we stop at this point the bucket-brigade approach would now 
determine the future behavior of the system by computing the matrix elements of the Hamiltonian 
between these two states and the metric formed by taking the overlaps of the two states. This 
result alone tells us that tunneling takes place. However, and that is what we now wish to study, it 
severely underestimates the tunneUng rate when the mass is large and the wells are well separated. 




(8.16) 



(8.17) 



(8.18) 



H = :^P^ + A(x2-/)^ 
2m 



(9.1) 
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The question is "what states do we have to add in order to compute the tunneling rate accurately 
?". The answer is, of course, determined by the instanton calculation. 

The key point is that in order to improve the energies of the two lowest states it suffices to 
apply the operator e to compute the effective Hamiltonian [|l|]-[^-[||] 

nij{t) = [{^ri\e-'^W,)Y''\^fi\e-i''ne-'^''W^) [{^f^\e-'''Wj)Y"\ (9.2) 

If we now follow the path integral procedure and divide the product up into a number of steps and 
insert a set of generalized Gaussian packets 

e-ivm+ip{t)^\y{t)) (9.3) 
where p{t) is chosen so that, as is always the case, 

pW = ifM. (9.4) 

and for simplicity y{t) is chosen to be a constant, then it follows that this contribution to the 
transition element is maximized if the function x{t) satisfies the equation 

^x{t)=AXx{t){x{tf-f) (9.5) 

where x{t) = — Cmin and dx{t = 0)/dt = 0, and 7 is chosen to be the same as that for the initial and 
final Gaussian packet. A picture of the solution and the discretized choice of a finite number of 
these states is shown in the appropriate section of the .html file. 

If we assume, as before, that these are the correct states to use to compute the time evolution 
of the initial Gaussian we simply compute the truncated Hamiltonian and exponentiate it. The 
corresponding animations in the accompanying .html file show that this computation is remarkably 
accurate for the real and imaginary part of the wavefunction, as a function of time, as well as for 
the tunneling rate. 

It is a straightforward matter to extend these ideas to field theory, however time and space 
preclude discussing this question at this time. 
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